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ABSTRACT 

By proposing a pure leptonic radiation model, we study the potential gamma-ray emissions from 
jets of the low-mass X-ray binaries. In this model, the relativistic electrons that are accelerated 
in the jets are responsible for radiative outputs. Nevertheless, dynamics of jets are dominated by 
the magnetic and proton-matter kinetic energies. The model involves all kinds of related radiative 
processes and considers the evolution of relativistic electrons along the jet by numerically solving 
the kinetic equation. Numerical results show that the spectral energy distributions can extend up to 
TeV bands, in which synchrotron radiation and s 5 mchrotron self-Compton scattering are dominant 
components. As an example, we apply the model to the low-mass X-ray binary GX 339-4. The 
results can not only reproduce the currently available observations from GX 339-4, but also predict 
detectable radiation at GeV and TeV bands by Fermi and CTA telescopes. The future observations 
with Fermi and CTA can be used to test our model, which could be employed to distinguish the origin 
of X-ray emissions. 

Subject headings: gamma rays: general — radiation mechanisms: non-thermal — stars: individual 
(GX 339-4) — X-ray: binaries 


1. INTRODUCTION 

Among Galactic X-ray binaries, there are about twenty 
strong candidate or firmly confirmed microquasars that 
present the extended relativistic radio jets, in which the 
jet is powered by an accretion of the central compact 
object. According to the mass of the companion star, they 
are classified as the high- and low-mass microquasars. 
For this kind of X-ray binaries, it is widely believed that 
during the low/hard spectral state, radio and infrared 
(IR) emissions originate from S 5 mchrotron processes of 
relativistic electrons in a persistent jet. 

In the low/hard state of these X-ray binaries, there 
is still ongoing debate regarding the origin o f X-ray 
emissions (e.g., IMcClintock & RemillardI 120061) . One 
possibility i s that they o riginate in a hot plasma 
corona (e.g., iLiu et al.l[l9^ . through the Comptoniza- 
tion of the thin accretion disk photons by hot elec- 
trons in the corona, or in a ho t accretion flow (e.g., 
lEsin et aT] 12001b IVuan et aO 120051) , via the Comptoniza- 
tion of S5mchrotron photons within the hot accre¬ 
tion flow. Another possibility is that they may 
have their orig ins from the jets, ei t her a s synchrotron 
emission (e.g., Markoff et al.l I2001L120031: iKaiseil 120061: 


Pe'er & C asellal ^091 IMait ra et all 120 091: 1 Vila & Romera 

20101 IPe'er & MarkoU 120121: IVila et all 120121) or s^ 


The study of observational and theoretical aspects of 
microquasars, particularly in the GeV and TeV domains, 
is a very active topic. From an observational point of 
view, five sources (including confirmed and candidate 
microquasars) are detected at high-energy (GeV) bands 
by the Fermi Large Area Telescope {Fermi LAT) and the 
AGILE satellite, and/or at very high-energy (TeV) bands 
by the Major Atmospheric Gamma Imaging Cherenkov 
Telescope (MAGIC), the High Energy Stereoscopic Sys¬ 
tem (HESS) and the Very Energetic Radiation Imag¬ 
ing Telescope Array System (VERITAS): Cy gnus X-1 
(e.g., [Albert et alJ 120071: iMalvshev et al.l 120131). Cvenus 
X-3 (e.g.. lAbdo et all I2nn9at ITavani et all 120091). SS433 
jBordas et al. 20151) . LS 5039 (e.g.. lAharonian et al.ll2005l 
Abdo et al. 2009b [). and LS I -t61° 303 (e.g.. [Albert et all 
20061 [Abdo et al.[l2009ij) . These gamma-ray detections 


chrotron self-Coinpton scattering (SSC, [Markoff et al.[ 
[20051 [Maitra et alJ[2009l). or external-rad iation-Compton 
scattering (ERG. [Pe'er & Markofj [20121) . Alternatively, 
the jet emitting di sk (e.g., Zhan g et alJ |2013[) and disk 
corona-jet models JOiao & Liu[[2015[) are also proposed 
to interpret X-ray emissions from black-hole X-ray bina¬ 
ries. 
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imply that extreme particle acceleration is at work in 
these systems. However, none of these detected gamma- 
ray sources is a low-mass microquasar. One possible 
reason is that these systems cannot exactly produce any 
GeV and TeV photons, as a result of the expected in¬ 
efficient effect of a low-mass companion star. Another 
possible reason is that one has not put forward an obser¬ 
vational plan to detect this class of object. However, they 
are emitting detectable gamma-ray signature. To our 
knowledge, almost no campaign detecting gamma-rays 
fro m a low-mass microq uasar is carried out (however, 
see [Bodaghee et al.[l2013|) . 

Based on leptonic and/or hadronic considerations, 
many models for jet emissions of high -mass micro¬ 
quasars are proposed in the lite r ature (e.g. .[Romero et alJ 
20031: [Bosch-Ramon et al.[[2006[: [Khangulvan et al.[[200^ 


Zdziarski et al.[l2014l: [Zhang et alJ[2014|) . In these mod- 

e8igenjjT?.)3tellar companion plays a key role, providing 
a seed photon field in a leptonic model or a target for 
proton-proton/proton-photon interactions in a hadronic 
one, even resulting in a cascade process. However, very 
little attention is focused on the low-mass microquasars. 
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in particular, for gamma-ray productions in a leptonic 
model. The ha dronic model for a low -mass microquasar 
is proposed in iRomero & Vilal (120081) . in which theoret¬ 
ical spectral energy distributions (SEDs) from cascade 
processes induced by relativistic protons, can well ex¬ 
tend into PeV bands. Subsequently, the further re¬ 
fined models (called as lepto-hadronic ones) are devel¬ 
oped to investigate the low-mass microquasars GX 339- 
4 (IVila &Romeroll20Tflh and XTE J1118+480 dVila et al.l 
I2012|) . In these lepto-hadronic models, radio and X-ray 
emissions are explained as an origin of the jet, in which 
SEDs could predict the detectable emissions at GeV and 
TeV bands by Fermi EAT and GTA. 

A recent search on neutrino em ission from a sample of 
six m icroquasars is presented in I Adrian-Martrnez et al.l 
ll2014|) by the ANTARES neutrino telescope. This study 
has provided upper limits of the possible neutrino 
fluxes, which constrains the scenario that the ratio of pro¬ 
ton to electron power is far less than 100 for a hadroni c 
jet model of low-mass microquasars dZhang et alj|2010l) . 
Besides, almost all the fits to GX 339-4 in the lepto- 
hadronic models ne ed the ratio of prot on to electron 
power close to 1 (e.g.. lVila & Romerol201(l|) . In a hadronic 
or lepto-hadronic model, protons have to be accelerated 
to about 10 PeV in order to excite hadronic cascade emis¬ 
sions. Such a high energy actually requires the existence 
of an extreme physical environment. Nevertheless, the 
requirement for electrons in a hadronic model is rela¬ 
tively loose, that is, the maximum ene rgy of electrons 
is cq r nmonly below Te V energies (e.g., iRomero & Vilal 
120081: IZhang et aT] 120101) . Howeyer, the simulations by 
particle-in-cell (PIC) methods demonstrate the presence 
of yery strong c oupling between electrons and shock- 
heate d protons dSpitkoyskyl 120081: iSironi & Spitkoysl^ 
12011|) . which implies that the dissipated energies are 
mostly conyerted into primary electrons. Purthermore, 
the recent hybrid PIC simulations indicate that it needs 
a stringent condition for an effectiye acceleration of pro¬ 
tons, that is, only when s hocks propagate parallel to th e 
upstream magnetic field dCaprioli & Spitkoyskyll2014|) . 

Giyen the aboye analysis and a lower radiatiye ef¬ 
ficiency of protons than that of electrons, the lep¬ 
tonic models appear to be strongly supported oyer the 
hadronic ones. Therefore, in this work we carry out 
a study about broadband emissions of low-mass mi¬ 
croquasars in the framework of a pure leptonic model. 
We want to know whether a pure leptonic model can 
also proyide the expected gamma-ray emissions as a 
hadronic one. It should be noted that there are some 
works on the jet emission o f X-ray binaries in the frame¬ 
work of the leptonic model (lKaiseij2006l:lPe'er & Casellal 
|200i. A common feature is that they considered the 
effects of radiatiye and ad iabatic coolin g of relatiyistic 
electrons. More concretely. iKaiseri (I2006|) focused on the 
flat synchrotron spectra of partially self-absorbed jets by 
considering the different jet geometries, such as ballistic 
and adiabatic scenarios. By using the electron spectral 
distribution with a Maxwellian shape at low energies 
plus a power-law tail at high energies, iPe'er & Ca.sellal 
ll2009^ studied the features of radiatiye spectra from the 
jet of X-ray binaries in detail and considered that a single 
acceleration episode is at work at the base of the jet. In 
these works, the models can produce multi-wayelength 
emissions ranging from radio to hard X-ray (or soft 


gamma-ray) bands and do not need to re-accelerate the 
electrons in the jet. 

In particular, except for including the radiatiye and 
adiabatic losses of relatiyistic electrons, we consider the 
eyolution of the relatiyistic electron along the jet by the 
kinetic equation. Purthermore, the dissipation region is 
extended to a large scale of the jet where the maximum 
energy of ultra-relatiyistic electrons can be constrained 
by the first-order Fermi acceleration mechanism. It is 
more obyious that the SEDs in our work can extend into 
the GeV and TeV ranges. As an example of application, 
we use our model to fit the low-mass microquasar GX 
339-4. The results show that our pure leptonic model 
can not only proyide possible explanations for the cur¬ 
rent obseryations, but also predict GeV and TeV band 
emissions. 

The paper is organized as follows. Descriptions for a 
leptonic jet model are presented in Section 2. Numerical 
results of the model are shown in Section 3. Section 4 is 
the fitting results to GX 339-4. Conclusions and discus¬ 
sion are giyen in Section 5. 


2. MODEL DESCRIPTIONS 

Like the case of actiye galactic nuclei, the mechanism 
of jet formation of X-ray binaries remains an open ques¬ 
tion. The seminal theoretical works suggested that large- 
scale magnetic fields anchored in an accretion disk or a 
black hole induce the generation of jets by means of 
a magneto-centrifugal mechan ism dBlandford & ZnajekI 
ll977t ' fMandford & Paynelll982^ . According to this main¬ 
stream thought, the jet in this work is considered as 
a Po 5 mting flux-dominated one in its irmermost region 
and then conyerts the magnetic energy into the matter 
kinetic energy at larger jet height. Howeyer, it is un¬ 
clear that for a Galactic jet, where and how it conyerts 
from the Poynting flux- to matter-dominated jets and 
how efficient the conyersion is. It is stated that MHD 
instabilities are possible trigge r mechanisms of the con¬ 
yersion ('e.g.. lSikora et al.l20n.5l) . These instabilities could 
produce the magnetic energy dissipation yia shocks or 
magnetic recormection, resulting in an acceleration of 
relatiyistic particles. 

In this study, we focus on the jet radiatiye proper¬ 
ties rather than its dynamical structure. We consider 
that only electrons are accelerated up to relatiyistic ener¬ 
gies and produce emissions. Howeyer, the cold-proton- 
matter and magnetic energies dominate dynamics of the 
jet. The geometry of a low-mass microquasar is the 
same as a high-mass one. Por the latter, the high-mass 
companion has significant effects in the aspect of radia¬ 
tiye output, eyen on dynamical structures of the system. 
Howeyer, in the current scenario the companion effect 
to the system is negligible. The geo metry of our jet radi - 
ation model is similar to Pigure 1 of IZhang et alJ (120141) . 
which was deyeloped to study the high-mass micro¬ 
quasar Cygnus X-1. 


2.1. Electron Distribution 

We assume that the steady-state electron distribution 
in a conical jet is goyerned by the following kinetic equa- 
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tion (e.g.. lZdziarski et'n]l2014tlZhang et al.ll2015l) 


dissipation region are considered as 


i||r|ftC2=N(y,z)l + d 


r|ftcN()<,2)L 


= Qin(y,z), 


( 1 ) 

where the first term represents spatial advection, corre¬ 
sponding to the divergence, V • vN, in a spherical coordi¬ 
nate, and the second term energy losses of electrons. Qm 
is the injection rate of relativistic electrons. N is the num¬ 
ber density of electrons per unit volume, as a function of 
the electron energy y and the jet height z from the central 
compact object. Lj is the bulk Lorentz factor of the jet. 



- l/Lj is the bulk velocity, and c is the speed 


of light. After introducing Nj,(y,z) = Fj/ijcdt'n(ztand)^N, 
andQin(y,z) = Qinn(ztan6)^TjPjcdt' (here, 0 is a half-cone 
angle of the jet). Equ ation (fill is changed into the fol¬ 
lowin g form (see also. lModerski et alJl2003l:lZhang et all 

IMl . 


dNy{y,z) 

dz 


dy 


dy 

N,,(r,2)^ 


Qm(y,z) 

c^jEj 


( 2 ) 


where Qm = dNy/dt'. The energy loss rate of relativistic 
electrons along the jet is given as 

cftTj^df'L 3z' 


where t' is the time in the jet co-moving frame. The 
second term of the right-hand side of Equation (|3) indi¬ 
cates the adiabatic loss of electrons. The total radiative 
loss rates of an electron, (dy/dt’)rad/ include synchrotron, 
its Comptonization, and ERC from the disk and corona 
photons. 

It is not clear that which mechanism is operat¬ 
ing to accelerate electrons up to relativistic ener¬ 
gies. It could be shock acceleration, s tochastic accel¬ 
eratio n or magnetic recormection (see iZdziarski et all 
12014 for some discussions), or shock interaction in a 
magnetic recormection site (iLazarian & Vishniad 119991: 

Ide Gouveia dal Pino & Lazarian 1 120051) . The diffusive 

shock acceleration involving the first-order Fermi pro¬ 
cess is often regarded as the most effective scenario 
in jets. However, this process requires that electrons 
are pre-heated up the value comparable to the energy 
of thermal ions. Eortimately, a quasi-Maxwellian elec¬ 
tron distribution with strong coupling with protons has 
been found in PIC simulations of collisionless relativistic 


dSpitkovskv 

120081 ISironi & Soitkovskyl |2011h and non- 

relativistic i 

the current 
leaf mildly 
electrons w 
small contri 
low-energy 
ergized, are 
jets of blacl 
On the oth 
ybr, which 
bution, are 

Riquelme & Spitkovskvll2'01l|) shocks. Eor 

study, jets are relativistic or should be at 
relativistic. In fact, the low-energy tail of 
ith a thermal distribution has a relatively 
bution to the SEDs. On one hand, a sharp 
cutoff, Vhr, from which electrons can be en- 
assumed in IZdziarski et al.l (12014) for the 
c-hole X-ray binaries in the low/hard state, 
ar hand, a hard injection spectrum below 
could approximate the Maxwellian distri- 
used to model a sample of blazars (e.g.. 

IGhisellmi & Tavecchiol 120091 IZhang et all 120121). Simi- 


larly to the latter, the relativistic electrons injected in the 


Qin(y, z) — Kq 


1 

y^yl7 + y‘^' 


(4) 


where ybr is the break energy of relativistic electrons, 
p and q are the energy spectral indices of the electrons 
below ybr and above ybr, respectively. The normalization 
constant of the electrons. Kg, is determined by 


Ke = 


^rePljetha 


nigC^l 


'ymin yfyiy+Yi 


ydy 


(5) 


where Lace = MaccC^ is the accretion power of the system 
via Roche lobe outflows of the companion star. Here, 
Mace is the mass accretion rate. The power of relativistic 
electrons, Lrei = 7]reikjet, accormts for a fraction of the jet 
power of Ljet = JljetLacc- We further consider qjet to be 0.1 
and 7]rei an adjustable parameter, ymin (assuming to be 
1) and ymax are minimum and maximum energies of the 
relativistic electron, respectively. 

We adopt a theory of shock acceleration to obtain elec¬ 
tron maximum energy, considering the balance between 
an acceleration rate and total loss rates. Eor a standard 
first-order Fermi acceleration mechanism, the accelera¬ 
tion rate of an electron in an ordered magnetic field B(z) 
is written as 

dy ^ 7]aceB(z)c 1 
dz mgC^ c/SjTj 


where r]ac is a parameter that characterizes an acceler¬ 
ation efficiency in a dissipation region and e is the ele¬ 
mentary charge. The magnetic fields evolve along the 
jet in the form of B(z) = BoZm/z. 


2.2. Radiative Processes 

Synchr otron emissions ar e calculated by using the for¬ 
mulae in ISpada et all (120011) . According to the radiative 
transfer equation, the resulting synchrotron luminosity 
in solid angle of 6Qj is given by 



l>i[.,{y,v')6Ny{y,z)dy 



(7) 


where 6/', is the pitch-angle averaged synchrotron emis- 
sivity and Zg is the optical depth of s 5 mchrotron self ab¬ 
sorption. 6Ny = Alj,6Qj/Oj is the number density of 
electrons in each solid angle. 

Eor numerical calculations of both SSC and ERC of 
disk and corona photons, we used Equation (14) of 
iZhang et all (I2014j) . In practice, we consider a spheri¬ 
cal corona with a radius of 3Rm surrounding the central 
compact object. Here, Rm ~ lORg is a truncation radius 
of the standard thin accretion disk, and Rg = GMbh/c^ 
is a gravitational radius of the black hole with a mass 
of Mbh- It is generally considered that the hard X-ray 
observations with a power-law spectrum is explained as 
Comptonization of accretion disk photons by hot elec¬ 
trons in the plasma corona. Eor simplification, the pho¬ 
ton field of the corona is assumed to be a power-law plus 
an exponential cutoff, Ug = Here, Eg is the 

maximum energy of the corona photon and the photon 
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index is 1.5 in this work. The normalization constant 
Kc are deduced by 


Kc = 


f E 


( 8 ) 


where Lcor = fc^-acc is the bolometric luminosity of the 
corona and /c is a parameter. 

We approxi mate the accretion disk as a planar geome¬ 
try following iSikora et all (|2013|1 . The energy density of 
disk photons in the jet co-moving frame at a height z is 
expressed as 


disk 


ST^GMaccMBH 

Anc 



(1 - j8jCOS0ex)^COS0ex 
(Z2 + R2)i^2 


dR, 

(9) 


where cos0ex = z/ Vz^ -t- R^. Rj^ and Rout are the in¬ 
ner and outer radii of the disk, respectively. The emis- 
sion spectra of an accretion disk are calculated following 
iFrank et al.l ll2002l) , together with an irradiation effect of 
the outer disk by itself irmer zone and/or a corona. De¬ 
tailed treatment of reprocessing of X-rays in the outer ac¬ 
cretion disk would involve many parameters and com¬ 
plex processes; we thus parameterize the reprocessed 
fraction of bolometri c flux with a factor of /out = 0.1 (see 
IGierlmski et al.ll2009l for more details). It is noticed that 
the temperature for a standard thin disk is T oc 
However, for an irradiation disk, it is T oc R Othe r 
model details are the same as those of iZhang et al.lll2014l) . 

In this study, we consider an extended dissipation re¬ 
gion that is located in the jet. In this case, an absorp¬ 
tion of high-energy photons can be avoided even if an 
emission region is very close to the cen tral compact ob¬ 
ject. As shown in lVila & Romerd (120101) for the low-mass 
microquasar GX 339-4, the absorption of TeV photons 
by infrared synchrotron photons is negligible even for 
a slightly compact dissipation region. From a research 
experience of high-mass microquasars, we know that 
at large jet scales absorptions of high-energy photons 
by disk and corona photons can be neglected. There¬ 
fore, the electromagnetic cascade processes due to the 
internal synchrotron photon, and the external disk and 
corona photons are not modeled. 


3. NUMERIGAL RESULTS 

In order to test our radiation model, we first explore 
the evolution of relativistic electrons and multiband 
SEDs by using some typical parameters for a general 
low-mass microquasar. For n umerical pro c edure s, in¬ 
terested readers are referred to lZhang et al.l (I2014|) . The 
used characteristic parameters for a general low-mass 
microquasar are given as follows: a mass of black hole 
of Mbh = 10 Mq, a distance of d = 2 kpc, and an accretion 
rate of Mace = 3.0 X 10“®Mq yr The model parameters 
are listed in Table[l]for Case A. 

The distributions of relativistic electrons as a function 
of the energy y along with the height of the jet z are plot¬ 
ted in Figure [l] From the beginning of the injection of 
relativistic electrons, the number of electrons increases 
with the height of the jet. It is noticed that when an 
injection spectrum of electrons is adopted in the form 
of a broken power law, the resulting electron spectra 
present the broken form around yb ■ At the low-energy 


limit, the number of electrons increases with their ener¬ 
gies. However, due to a rapidly cooling of electrons, the 
number of electrons decreases with their energies at the 
high-enemy limit. 

Figure |2] demonstrates broadband SEDs for a general 
low-mass microquasar. As shown in this figure, the 
emissions from the disk (including a standard thin disk 
and the corresponding irradiated disk) and hot plasma 
corona contribute radiative fluxes ranging from near-IR 
through ultraviolet (UV) to MeV bands. The synchrotron 
emission provides radiative fluxes at radio, IR, and about 
0.1 GeV bands. The SSC spectrum dominates fluxes 
above 1 GeV. Even though relativistic electrons injected 
in the jet is at a distance of 70Rg close to the central 
compact object, the ERG of soft photons of the accretion 
disk and corona are still negligible. 

As a result, we find that the nonthermal emissions 
from the jet, e.g., synchrotron and SSC, dominate the 
jet radiative outputs. It is well-known that these two 
components are mainly governed by electron distribu¬ 
tions and a magnetic field strength in an emission re¬ 
gion. Therefore we define a magnetization parameter 
a = U’e,lu\an to characterize a role that magnetic fields 
play in a jet flow, that is, a ratio of magnetic energy den¬ 
sity to matter kinetic density. Here, mb = B^/Sn is the 
magnetic energy density in a completely ordered mag¬ 
netic field and Wbin = Ljet/cFj/3j7T(ztan0)^ is the kinetic 
energy density of jet matter. With the values adopted 
in this numerical test, we find cr ~ 4 x 10“^ at the initial 
location of the dissipation region. It is very evident that 
the entire dissipation region is sub-equipartition. 

4. APPLICATION TO GX 339-4 

GX 339-4 has been intensively obse rved in ra- 
dio, IR, optical, and X-ray bands (e.g., ICorbel et alj 
200(il: IHoman et all 1200.51: lYu et all l20(M IWu et all 1201(1 
Gandhi et all 120111 Tv^& Yul 12012ft . The source goes 
through all the spectral states of a typical black-hole 
X-ray binary. Tbe radio image of the GX 339-4 jet 
are achieved by IGallo et'aLl (120041) after the X-ray out¬ 
bursts of 2002, which confirmed that this source is a mi¬ 
croquasar. A central black hole, with a mass > 6Mq, 
is in an orbit aroimd a dark com panion star, with a 
mass in the range [0.166Mo,1.1Mq] (iMunoz-Darias et al.l 
120081) ■ which is loca ted at a distance of ~ 6-9 kpc 
dZdziarski et alj|2004ft . 

It can be seen that this source is a low-mass micro¬ 
quasar. A distance of d = 6 kpc and a black-hole mass 
of Mbh = 6 are used in our following fittings. In addi¬ 
tion, some typical values are fixed as a jet bulk Lorentz 
factor of Fj = 2, an accretion rate of Mace = 5 X 10“®Mq 
yr“^, a half-cone angle of 0 = 5°, and a viewing angle 
of 0obs = 30°- As exposed in Sectionj^ s 5 mchrotron and 
SSC dominate the radiative outputs of the jet for the case 
of the sub-equipartition. Below, we are to explore two 
cases, that is, the magnetization parameter cr ~ 1 and 
<< 1 . 

The first case corresponds to an equipartition between 
magnetic energy density and kinetic energy density of 
matter. The initial location of the emission region is first 
fixed at a distance of 10® cm (~ 67.3Rg). We then ad¬ 
just the values of the jet power Ljet and the magnetic 
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Fig. 1. — Distributions of relativistic electrons as a function of electron energy at different heights of the jet. The adopted parameters are listed in 
TableUfor Case A. 



Fig. 2.— The SEDs of a general low-mass microquasar. The adopted parameters are listed in Table[T]for Case A. The disk spectrum includes the 
contribution from an irradiated outer disk (short dashed line). 
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Fig. 3.— Fitting the SEDs of GX 339^ for the c orona-dominated X-ray emissions under an equipartition condition. The radio through UV up to 
X-ray observations are from IGandhi et all 12011]) . The thick solid lines are 5a sensitivity limits for different instruments (1 yr survey mode for Fermi 
and 50 hr of direct exposure for HESS, MAGIC-II and CTA). The fitting parameters are listed in Table[T]for Case Bl. 
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Fig. 4.— Fitting the SEDs of GX 339^ for the corona-dominated X-ray emissions under a sub-equipartition scenario. The thick solid lines are the 
sensitivity limits of Fermi, HESS, MAGIC-II and CTA. The fitting parameters are listed in Table[l]for Case B2. 
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field strength Bq to make a close to 1. At last but no 
least, some parameters that are associated with relativis¬ 
tic electrons, such as spectral indices and break energy, 
are adjusted to fit observations. According to the con¬ 
straints of X-ray observations, which are regarded as 
the corona-dominated emissions, the value of the break 
energy ybr is determined. 

Figure|l]shows the SEDs of GX 339-4 confronting with 
the observations available and the sensitivity limits of 
telescopes. The fitting parameters are listed in Table[T]for 
Case B1. As shown in this figure, the synchrotron emis¬ 
sion and the irradiation of the accretion disk can well 
explain observations ranging from radio to UV bands. 
X-ray data are fitted by the corona photon spectrum that 
is considered as a power-law plus an exponential cutoff. 
Synchrotron and SSC spectra are dominant at the GeV 
and TeV bands, respectively. It can be seen that syn¬ 
chrotron emission flux can compare to the sensitivity 
limit of Fermi, however, SSC flux is significantly below 
the sensitivity limit of CTA. The ERG spectra from disk 
and corona are negligible. 

By adopting the similar method mentioned previously, 
a sub-equipartition case is explored here. The fitting re¬ 
sult is presented in Figure |4] with the parameters given 
in Table [T] for Case B2. From this figure, we see that at 
low energy bands, the model can well reproduce obser¬ 
vations. SSC component dominates the radiative flux at 
GeV and TeV bands; this flux can match the sensitivity 
limits of Fermi and CTA, but still below the sensitivity 
limits of MAGIC-II and HESS. In this case, we obtain a 
magnetization parameter a = 3.44 X 10“® at the initial 
location of electron injection, which indicates that the 
emission location is in a very low magnetization (e.g., 
matter-dominated) region. 

Hereinbefore, we consider that the X-ray emission has 
its origin in a corona. However, studies of correlation 
between X-ray and radio flux show the fact that there 
is a tight positive correlation for GX 339-4, Fradio F°-^ 
jCorbel et al.ll2003h . which appears to suggest that both 
waveband emissions have a common origin, that is, they 
may be from the jet region of X-ray binaries by possible 
radiative processes (iCorbel & Fendedj 2002l:ICorbel et all 
2nnl IMarkoff et a 111200111200.31.1200.31. iMaitra et al JI2n(m 
Vila & Romeroll20Tc , see also Section 1 for a review). Af¬ 
ter weakening the contribution from the hot corona, we 
use the jet nonthermal emission to fit X-ray observations 
together with other band data, considering two scenar¬ 
ios, e.g., equipartition and sub-equipartition. 

In order to fit X-ray observations, we have to adjust 
the spectral index q below 1.5. Meanwhile, we must 
use a more low break energy ytr close to the minimum 
energy (set as 1) in order to fit IR and optical data. In 
view of this reason, we give up using the broken power- 
law spectrum and assume a single power-law form, 
Qin oc y P, ranging the entire energy band. Figures |3 
and [6] show the fitting SEDs of GX 339-4 for equiparti¬ 
tion and sub-equipartition cases, respectively. The used 
parameters are listed in Table [T] for Case Cl (Figure |5) 
and Case C2 (Figure |6]l. As shown in these figures, the 
synchrotron emission can reproduce radiative fluxes at 
radio, IR, optical, and X-ray bands. SSC component 
dominates gamma-ray band emissions, which can pre¬ 
dict detectable signature by Fermi and CTA telescopes. 


but still below the sensitivity limits of MAGIC-II and 
HESS. It is noticed that for the case of equipartition the 
synchrotron emission presents a higher peak frequency 
than that of the sub-equipartition one (a = 2.39 X 10“^), 
because a stronger magnetic field is used in the case of 
the former. 

In summary, we have employed our leptonic radia¬ 
tion model to fit observations from GX 339-4 and predict 
possible gamma-ray emissions. For the two equiparti¬ 
tion scenarios (shown in Figures |3] and |5]l, it needs to 
provide a high acceleration efficiency of J]ac = 0.1. How¬ 
ever, this value is relatively low, that is, 0.01, for the 
sub-equipartition cases (see also TableH). The main rea¬ 
son is that the expected maximum energy of relativistic 
electrons is required close to about 1 TeV in order to emit 
gamma-rays. In addition, because the dominant (syn¬ 
chrotron) loss rate is oc y^B^ and the acceleration rate 
oc r]acB (see Equation[6), we have the electron maximum 
energy of pmax (7?ac/S)^^^. The fittings show that Fermi 
and CTA should have abilities to detect this source. We 
note that synchrotron and SSC spectra present different 
slopes for the different origins of X-ray emissions. Thus, 
the possible detection with Fermi and CTA could help to 
distinguish the origin of X-ray emissions. 

5. CONCLUSIONS AND DISCUSSION 

We have proposed a pure leptonic radiation model 
to study electromagnetic radiation spectra from jets of 
the low-mass X-ray binaries. This model considers the 
evolution of relativistic electrons along the jet by numer¬ 
ically solving the kinetic equation. The main photon 
fields for a general low-mass microquasar, such as ac¬ 
cretion disk and corona photons, are involved in the 
model. Our results show that both synchrotron and 
SSC spectra play a dominant role, which is similar to 
the synchrotron plus SSC model that is often used to 
interpret the emissions of extragalactic BL Lacertae ob¬ 
jects. The ERC spectra (disk and corona) are negligible, 
though the injection location of electrons is close to the 
central compact object. As an example, we have applied 
this model to the low-mass microquasar GX 339-4. Our 
model can well reproduce the observations available at 
radio, IR, optical, UV, and X-ray bands, and predict po¬ 
tential gamma-ray emissions at the GeV and TeV bands. 
Observations with instruments, Fermi and CTA, can be 
used to test our model, which could be adopted to dis¬ 
tinguish the origin of X- ray emissions. 

Based on the work of IMarkoff et al.l II2005I) . an irradi¬ 
ated disk plus jet model was used to fit observations 
of the low-mass X-ray binaries, XTE T1118-h480 and GX 
339-4, from radio through X-rays ijMaitra et al.l 1200^ . 
Their results on the fitting of GX 339-4 showed that 
the jet model emitting synchrotron and SSC radiation 
dominates emission output throughout the entire wave¬ 
band, without requiring a disk component, in which SSC 
component provides a significant fraction of the X-ray 
emission flux especially when the base of the dissipa¬ 
tion region is set to close to the central compact object. 
In the current work, SSC spectrum dominates the GeV 
and TeV band emissions and has no contribution to X- 
ray observations. In addition, an irradiation disk and 
jet synchrotron or corona components are alternatively 
adopted to explain emissions from radio to X-ray bands. 
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Fig. 5. — Fitting the SEDs of GX 339^ for the jet-dominated X-ray emissions under an equipartition condition. The thick lines are the sensitivity 
limits of Fermi, HESS, MAGIC-II and CTA. The fitting parameters are listed in Table[T]for Case Cl. 



Eig. 6.— Fitting the SEDs of GX 339^ for the jet-dominated X-ray emissions under a sub-equipartition scenario. The thick lines are the sensitivity 
limits of Fermi, HESS, MAGIC-II and CTA. The fitting parameters are listed in Table[T]for Case C2. 
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As shown in Figures 3-6, this work uses more complete 
IR, Optical and UV data compared with that of Maitra et 
al. (2009). However, SSC component is still loose due to 
lacking any data above the X-ray band while the current 
observations can constrain the S5mchrotron spectrum. 

For the case of the corona-dominated X-ray emissions, 
locating the dissipation region at sub-equipartition 
scales is optimistic to gamma-ray productions. More¬ 
over, the value of pbr ~ 10^, which is needed to elimi¬ 
nate the synchrotron spectral contribution to X-ray emis¬ 
sions, is close to the average energy of the shocked- 
ions. For the case of the jet-dominated X-ray emissions, 
we used an electron spectrum with a single power-law 
form. However, in this case, a hard electron spectral 
index of p = 1.4 is used. In the theory of non-relativistic 
diffusive shock acceleration, it is well-known that the 
spectral index should be close to 2. At present, though 
there are many works on this topic, we still know little 
for more physical aspects. We notice that a particle in- 
iec tion spectrum with a hard index has been obtained 
in IStecker et al.l II2007I) . by using a Monte Carlo sim¬ 
ulation of particle diffusive acceleration at relativistic 
shocks. Besides, formation of hard power-law distribu¬ 
tion 1) of the energetic particle spectra is derived by 
3D P IC simulations f rom relativistic magnetic reconnec¬ 
tion jGuo et al.ll2014D . 

From the perspective of energy spectral fitting, many 
works regarding the jet emissions indeed demon- 
strat e that the power law index is greater than 2 

te.g..lMarkoff et ^1^051: fKaisedl^nO^lMaitra et alJl2009l: 
iPe er & Casellal 120091) ■ in which almost no gamma-ray, 
especially very high-energy gamma-ray emission is ex¬ 
pected. However, if one expects a pure jet model (espe¬ 
cially in a low-mass X-ray hinary) to produce the entire 
broadband emissions ranging from radio via optical and 
X-ray to GeV and TeV bands, the jet itself needs to pro¬ 
vide more energy to the energetic electrons. In this case, 
a hard power-la w index is needed, as shown in Figures 
5 and 6 (see also lVila & Rome^l2010l:lZhang et al.ll20l4 
for a similar scenario), though one would expect a large 
power law. 
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As pointed out in Section |2l an expected location of 
particle acceleration should not be at the base of the jet, 
in which outflow should be a Poynting flux-dominated 
one. According to the fitting results of GX 339-4, it ap¬ 
pears also to support the scenario that the jet emission 
is away from the central compact object, and at a sub- 
equipartition region (see also Figures |4] and [6]l. In the 
present work, we parameterized the acceleration effi¬ 
ciency as a factor rjac- Without a doubt, this overlooks 
many physical details regarding jets and precludes in¬ 
sight into the jet dynamical properties. Furthermore, 
in order to obtain the expected maximum energy of 
relativistic electrons, which mainly decides the condi¬ 
tion whether electromagnetic spectra can extend into 
gamma-ray regimes, we have used the relatively effec¬ 
tive efficiency ijac = 0.1 or 0.01. 

For the high X-ray luminosities observed in black-hole 
X-ray binaries or extragalactic blazars, if their emission 
fluxes are produced in a jet, a powerful jet with a sig¬ 
nificant amount of leptonic content is necessary, because 
electrons have very high radiative efficiency in contrast 
to protons. A lepto-hadronic model needs more effective 
acceleration efficiency than a pure leptonic scenario, that 
is, it is necessary to make protons accelerate up to about 
10 PeV, at which protons can efficiently excite proton- 
photon/proton-proton interactions. In practice, because 
the lepto-hadronic model involves multiple secondary 
particles promoting cascade processes, it is difficult to 
distinguish the contribution of each individual compo¬ 
nent. In this regard, our leptonic radiation model is more 
pure. 
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TABLE 1 

The Free Parameters Adopted in the Model. 


Case 

Biel 

^m[^g] 

^end [^g] 

Bo[G] 

Tbr 

P 


^ac 

/c 

£c[keV] 

A 

0.1 

7.00 X 10 ^ 

3.00 X lO *” 

1.0 X 10 ® 

1.0 X 10 ® 

1.0 

2.5 

0.01 

0.01 

100 

B1 

0.009 

6.73 X 10 ^ 

3.37 X 10 ^ 

2.0 X 10 ® 

2.0 X 10 ® 

1.0 

2.6 

0.1 

0.12 

60 

B2 

0.08 

6.73 X 10 ® 

4.04 X 10 ® 

3.0 X 10 ® 

1.0 X 10 ® 

1.0 

2.6 

0.01 

0.12 

60 

Cl 

0.07 

6.73 X 10 ^ 

3.37 X 10 ^ 

3.0 X 10 ® 

- 

1.4 

- 

0.1 

0.001 

60 

C2 

0.04 

6.73 X 10 ^ 

4.04 X 10 ® 

2.5 X 10 ^ 

- 

1.4 

- 

0.01 

0.001 

60 


* Note. Symbol indicating i/reT conversion efficiency of relativistic electrons; Zin: start of dissipation 
region in jet; z^nd^ end of dissipation region in jet; Rg: gravitational radius; Bg: magnetic field strength; 
ybr: break energy of electron; p: spectral index of electron below spectral index of electron 

above ybrl ^ac: acceleration efficiency; fy. ratio Lcor/Lacci Sc'- break energy of corona photons. 
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